Fixed Points of the Dissipative Hofstadter Model 
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The phase diagram of a dissipative particle in a periodic potential and a magnetic field is studied 
in the weak barrier limit and in the tight-biding regime. For the case of half flux per plaquette, and 
for a wide range of values of the dissipation, the physics of the model is determined by a non trivial 
fixed point. A combination of exact and variational results is used to characterize this fixed point. 
Finally, it is also argued that there is an intermediate energy scale that separates the weak coupling 
physics from the tight-binding solution. 
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Introduction A quantum particle interacting with an 
environment with a macroscopic number of degrees of 
freedom, the Caldeira-Leggett (CL) model Q, is one of 
the simplest models used in the study of decoherence in 
quantum systems. This model has been generalized to 
include the motion of a dissipative particle in a periodic 
potential Q , in a finite magnetic field Q , and in a com- 
bination of both situations |^y. The problem described 
by such a model applies to a large number of situations 
in condensed matter, quantum computation, and string 
theory. A few examples are: flux qubit dephasing in 
quantum computers [3, defects in Luttinger liquids 0, 
junctions between many Luttinger liquids 7|, and non- 
trivial backgrounds in open string theory ]4$. The CL 
model is also relevant to the study of the dephasing in- 
duced in mesoscopic systems by external gates |g, |9j , and 
it always reproduces the short time dynamics of particles 
interacting with ohmic environments |lCj . The model 
also describes the quantum motion of a vortex in a lat- 
tice. This model has attracted interest in the study of 
d-wave superconductors with strong phase fluctuations 
[ilf . Note that, in this context, dissipation due to low 
energy modes arises naturally. 

Although the phase diagram of a dissi pati ve p article in 
a periodic potential is well understood [l^, UM , there is 
no similar degree of understanding when a magnetic field 
is also added to the problem0,Q|. This model became 
known as the dissipative Hofstadter model, and in the 
present work we analyze its phase diagram with a square 
lattice symmetry, the renormalization group (RG) flows 
and fixed points. We use mappings into spin and fermion 
Hamiltonians, and variational methods in order to ob- 
tain further information on the phase diagram, which, as 
discussed below, presents a number a new features with 
respect to the model without a magnetic field. 

The model. In the absence of dissipation, the Hofs- 
tadter problem has a complex energy spectrum |lfj| . The 
model has a duality between the weak coupling and tight 
binding limits of the periodic potential and, in both 
cases, the spectrum can be described by Harper's equa- 



tion |ljj • A similar duality holds in the dissipative case 
0, which admits further extensions (see below). 

We start by considering the limit where the periodic 
potential is weak. It was shown in Ref. that from 
perturbation theory on the lowest Landau's levels the 
dissipative model can be described by a boundary con- 
formal field theory in (1 + 1) dimensions with action (we 
use units such that h = 1 = &_b), 
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where is the totally anti-symmetric tensor. The parti- 
cle's coordinates are represented by the boundary degrees 
of freedom of the field, R(r) = 6i(0, r)e x + 62(0, r)e y . 
In addition, the particle moves in a periodic potential of 
lattice spacing a and amplitude A (A = V/A -C 1, where 
V is the potential strength, and A is a high energy cut- 
off), and is subject to a perpendicular magnetic field of 
amplitude j3 = Ba 2 /&o (where $0 is flux quantum). Dis- 
sipation arises from the first term in the r.h.s. of Eq. JJJ 
when the bulk modes (Oj(x, r) with x > 0) are traced 
out. The dissipation strength is given by a — rja 2 (where 
rj is the dissipation coefficient). 

In the absence of the potential the theory is Gaussian 
and the field propagator reads 0, 0, Q] 



(e i (0,r)6 J (0,0)) = 2aU\r\Si, 



iir<3 sgn (r) e zJ ,(2) 



where a = a/{a 2 + (3 2 ) and (3 = f3/{a 2 + (3 2 ). The 
first term in Eq. J5J) is the well studied logarithmic cor- 
relations. The second part of the propagator is the 
Aharonov-Bohm phase that the particle picks due to the 
magnetic field. Using this result, it is simple to write the 
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partition function as an expansion in powers of A 

z = E E E A " 
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where A ± (r,) = e ±4e i(°> T ;). Eq, (J^J has a simple physical 
interpretation: each insertion of represents a jump of 
the center of the particle's Landau orbit by a vector of the 

; +na/ \J a 2 + f3 2 e y 
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"dual" lattice f mi „ = ma/ 

(n and m are integers) |l4| . Which, for a = is a distance 
proportional to the Larmor radius (w c = /3 _1 / 2 <z). 

The complementary limit to the physics of Eq. is 
to consider very large barriers. Thus, instead of the low- 
est Landau orbits, a tight binding approximation to the 
spectrum in the absence of dissipation is natural start- 
ing point0,Q. The partition function is now expanded 
in powers of the nearest neighbor hopping amplitude, t, 
between the minima of the periodic potential. The result 
is identical to Eq. © with the substitutions of Table |U 
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tight binding 
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particle in a triangular lattice can have a non trivial fixed 
point at intermediate coupling. 

Variational treatment. For g < 1, f2(£) scales to- 
ward strong coupling. This usually suggests that the 
fields x ,y(0) become "pinned" at some value 9 XtV (Q). We 
can gain insight into this "pinned phase" usin g th e Self- 
Consistent Harmonic Approximation (SCHA) 20]. This 
approximation replaces the original periodic potential by 
harmonic wells adjusted self-consistently. Within SCHA 
we replace the boundary term in Eq. (2J by 



E 



K a [6a(P) 



where K x and K y are variational parameters. The p- 
dimensional state |0) has also to be adjusted variation- 
ally, and 6 x ,y = arg(0|T X)y |0). 




TABLE I: Duality relations between the strong and weak 
coupling limits of the dissipative Hofstadter model. 

In the following, we focus on a cross section of the 
phase diagram, /3 or j3 = 1/q with q 6 Z, that contains 
most of the interesting features. We found convenient to 
re- write the problem in an unified Hamiltonian formalism 

n = | jH ' dz{j-[d z e x Jz)] 2 + 2g[Tl x ^z)] 2 } 
+ AQT^^ + AnT/'M+Lc, (4) 

where we set v = 1, [9 x>y (zi) , U Xty (z 2 )] = 8 x , y 5 {z\ - z 2 ) 
and T Xj y are p-dimensional matrices (%^y = Tj y ) that 
satisfy the algebra 

r x r y = e 2 -/? Ty r x . (5) 

The correspondence between the parameters of Eq. l|4"| 
and the dissipative model are summarize on TableQ] The 
stability of the both limits is given by the lowest order 
renormalization group (RG) equation 

d e n = (i - g )n , (6) 

where d£ = dA/A. Since the scaling dimension of fl(£) 
in the strong coupling case is not the inverse of the one 
at weak coupling, there are values of (a, (3) where both 
A and t have runaway flows. This is similar to the case 
considered in Ref. 0, where it was shown that a 



FIG. 1: Set of degenerate minima obtained using the SCHA 
for q = 2. The arrows denote the orientation of the state 
\m,n) for m,n = 1,2. Possible hopping terms between these 
minima are also sketched. 

Because of the periodicity of the potential, states |n) 
such that (n|T Xj y|n) differ by a phase, have the same 
ground state energy. It can be shown that the lowest 
energy is obtained for |(0|T X |0)| = |(0|7y|0)| = T and 

K x = K y = A QT (QT)^ 1 . (7) 

Given a state |0), and using the relations in Eq. JSJ, 
we can construct p 2 states, \m, n) = T x n Ty\0),m,n — 
1, • • • ,p — 1 that lead to the same energy. Thus, the 
SCHA leads to a degenerate set of states labeled by 
the minimum in the periodic potential. The situation 
is illustrated in Fig.p^, where the case q = 2 is shown. 
The scaling dimension of the hopping between minima 
in the sublattice defined by a given vector \m, n) is 1/g. 
Although tunneling is also possible between minima in 
different sublattices, it is reduced by the overlap factor 
| (m, n\m! , n') |. These minima are closer in real space, and 
their scaling dimension is 1/ (Ag). Hence, for 1/4 < g < 1 
both the weak and strong coupling limits are unstable 
further supporting the existence of an intermediate fixed 
point. For a general value of (3 (or (3) — 1/p this result 
generalize to 1/p 2 < g < 1. Moreover, the overlap of 
the states |m,n) define an effective Berry's phase gen- 
erated when moving around each plaquette. It is easy 
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to show that the flux per plaquette needed to generate 
this Berry's phase is p. Hence, within the SCHA the 
weak barrier problem has an additional duality property, 
which leaves (3 ((3) unchanged, but replaces g <-» l/(p 2 g). 

Mapping to a spin chain for q = 2. For [3 = 1/2 in 
the weak coupling case, or = 1/2 in the tight binding 
limit, the operators T x and T y in Eq. J3J reduce to Pauli 
matrices, a x ,a y . This equivalence suggests the use of a 
spin Hamiltonian with the same universal properties for 
the environment. Thus we replace Eq. by two semi- 
infinite XXZ chains 

ft = E °n°n+l + CT n°"ra+l + ^nK+l 

+ v s n{^o+^i) (8) 

where A = cos[7r(l — g)] and v s — ir\ sm(7rg)|/(7r + 
2 arcsin[cos(7rg)]). 

The spin Hamiltonian provides a different perspec- 
tive of the infinite coupling limit studied by SCHA. AS 
O — > oo, the low energy sector tends towards the tensor 
product of two semi-infinite chains (starting from sites 
±2) plus the low energy excitations of the three strongly 
coupled spins at sites -1, and 1. As the SCHA sug- 
gested, there are four degenerate states (see Fig. (Q). 
When we consider < oo, the interaction between sites 
±2 and ±1 can be treated as a perturbation of order 
I/O and we the degeneracy is lifted to a doublet. In fact, 
this doublet is protected by a hidden symmetry (see be- 
low). After defining dual spin variables [2l|, = cr^er* +1 
and t* — Y\j <n Vj, we find that [tq,7Y] = 0. This con- 
served quantity is non-local in the original spin (dissipa- 
tive) problem, thus it correspond to a topological charge. 

We can further understand the intermediate fixed 
point by solving the "non-interacting" problem, g — 
1/2. As a bonus to be solvable, it is also believed that 
this point separates four different phases in the (a, (3) 
plane0. Using the dual spin variables, the Hamiltonian 
breaks into three independent parts, 

ft = E i< + <-i<] + E [a£ + + V , (9) 

n^O n 

with the definitions: \i\ = r\ n , v l n — r\ n _ x and V = 
(v s rt— 1) [fJ.Q + M-iMo] • Eq. JHJl implies that the odd sites 
of the original chain are mapped into two semi-infinite 
quantum Ising chains with open boundary conditions. 
The even sites are mapped into a single quantum chain 
and an impurity term (V). After fermionizing the three 
chains and taking the continuous limit, it is straightfor- 
ward to show that V is an irrelevant operator of dimen- 
sion 2. Hence, the manifestly conformal invariant RG 
fixed point is v s fi = 1. In the fcrmionic language, the 
conservation of the topological charge is represented by 
a single Majorana fermion localized at the origin. In ad- 
dition, we just showed that for g — 1/2 the fixed point 
is the resonance condition to a fermionic channel. Since 



dim V = 2 at the "non- interacting" point, it is very likely 
that V will also be an irrelevant operator for other val- 
ues g. Eventually, as we consider g — > 1, the repulsive 
interaction between fermions became sufficiently strong 
to close the fermionic channel through the localization of 
a second Majorana (O — > fixed point). 

The correspondence with the SCHA give us as simple 
picture about the particles mobility. In SCHA, the four 
minima of the potential are organized in sub-lattices de- 
picted in Fig. JIJ. With the mapping to the spin chain, 
they can also be classified accordantly to the two possible 
values of the topological charge (t§ ) . This fact suggests 
that at the intermediated fixed point the lattice breaks 
into two sub-lattices. Tunneling between minima of dif- 
ferent sub-lattices does not occur, while the amplitude 
for wells in the same sub-lattice is given by the renor- 
malized value of I/O. This is very similar to the inter- 
mediated fixed point of a Brownian motion in a trian- 
gular lattice^^j where there are three geometrical sub- 
lattices. For intermediate values of dissipation, there is a 
regime where the particle avoids one of the sub-lattices, 
but moves on the other two. This scenario of an interme- 
diated mobility can be further supported by noticing that 
for the exact solution, g = 1/2, the current operator also 
becomes quadratic in the fermion operators. Since the 
correlations will decay as r~ 2 at long times, the particle 
mobility, pij = lim^^o ^(©i(0, cj)8j(0, — u))), is finite at 
the fixed point |2^ . 

Fixed point at or near g = 1. When a, a = 1 the 
diagonal correlations between the e 10 *^ , e l8v ^ operators 
decay as r -2 . The RG equation can be derived in an 
e= 1-50123 

expansion scheme, 

d e n = ett-C sin 2 (w/q) O 3 + O (e 2 , O 5 ) , (10) 

where C is a constant of order unity. For q = 2, Eq. (|10|) 
implies a renormalized f2 cx \fl. As g — > 1, this fixed 
point merge with the trivial O = 0. The physical mean- 
ing of O = is straightforward when we look from the 
perspective of a quantum impurity problem. In Eq. Q), 
the two bosonic fields favor the localization of the spin 
variable along orthogonal directions. Thus, when g = 1, 
the "frustration" decouples the spin from the baths|24|. 

Phase diagram for half flux per unit cell. We now focus 
on a magnetic filed which corresponds to half flux per 
plaquette, (3 = 1/2, which illustrates the different fixed 
points mentioned above. We summarize the discussion 
on Fig. J2J). 

In the tight binding limit, the particle is localized for 
a > 1. Close to a £ 1 there is an intermediate fixed 
point, t* cx s/l — a. For a ~ 1/2 the exact solution 
shows that t* ~ 1. Finally, the duality transformation 
obtained by variational means indicates that a localized 
solution is unstable for a > 1/4, so that an intermediate 
fixed point exists for 1/4 < a < 1. 

In the weak barrier limit, except at a — 1/2, the A = 
fixed point is unstable for all values of a. For a = 1/2 the 
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(a) weak barrier limit. (b) tight binding limit. 

FIG. 2: Phase diagram of the dissipative Hofstadter model 
with half flux per unit cell, /3 — 1/2. The lines and symbols 
show the expected fixed points. 

RG flow of A is zero, leading to a line of fixed points . 
This point in the phase diagram is equivalent to the well 
known line of fixed points of the model without the mag- 
netic field@|. This happens because as the particle hops 
in the "dual" lattice, r, it picks a phase of 2n around each 
plaquette. For a = y/3/2, the partition function is iden- 
tical to the partition function in the tight-biding limit. 
Since the model is self-dual there ought to be at least one 
fixed point at intermediate coupling. For a = l/(2\/3) 
we find f3 = 3/2. The model has the same properties as 
when j3 — 1/2, and the action, eq.JD, is formally equiva- 
lent to the action obtained in the tight binding limit. It 
seems likely that the fixed point obtained from the vari- 
ational approach in this regime has the same properties 
as the one in the tight binding limit. 

In the region < a < 1/4 the weak and tight binding 
limits have RG flow towards strong coupling. The SCHA 
suggests that the particle is indeed delocalized with a 
phenomenology quite different from the A = fixed point. 
Instead of moving in the "dual" lattice, r, it freely moves 
in the lattice induced by the potential (23 . 

The existence of the self-dual point strongly sug- 
gests that for some parameters both the weak and the 
tight binding limits can be used to describe the model. 
However, at a = 1/2 the different approaches lead to 
markedly different results. Similar discrepancies do exit 
in many other parts of the phase diagram. Hence, it is not 
obvious how to extrapolate the results from the weak bar- 
rier case to the tight binding limit and vice-versa. These 
differences between the weak barrier and tight binding 
limits are related to the range of validity of the field the- 
ories that describe each one. For instance, when the RG 
flow of the weak coupling case leads to energy scales of 
the order of max(?7,a; c ), Eq. @J is no longer justified. 
Then, the theory must be supplemented with operators 
due to transitions to higher Landau levels. This is clear 
in the a = 1/2 case, where the particle in the weak bar- 
rier limit effectively hops in a "dual" lattice with lattice 



parameter \[2a. Hence, starting from Eq. it is not 
possible to account for the effects of the particle tunnel- 
ing between minima of the periodic potential separated 
by a. This means that at a certain energy scale the line of 
fixed points stops, and the problem starts to renormalize 
to the exact solution that we discussed in the text. Con- 
versely, using the duality properties of the model, there 
are other regions of the phase diagram where the tight 
binding suffers by the same problem. 

In conclusion, we studied the dissipative Hofstadter 
model using scaling, exact results, and a variational ap- 
proach. This allowed us to characterize the intermediate 
coupling fixed point of the model. Finally, we showed 
that results obtained in weak barrier or tight binding 
limits cannot be straightforward connected. 
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